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We present a mathematical procedure which leads us to obtain analytical solutions for the atomic 
inversion and Wigner function in the framework of the Jaynes-Cummings model with an external 
quantum field, for any kinds of cavity and driving fields. Such solutions are expressed in the integral 
form, with their integrands having a commom term that describes the product of the Glauber- 
Sudarshan quasiprobability distribution functions for each field, and a kernel responsible for the 
entanglement. Considering two specific initial states of the tripartite system, the formalism is then 
applied to calculate the atomic inversion and Wigner function where, in particular, we show how the 
detuning and amplitude of the driving field modify the entanglement. In addition, we also obtain 
the correct quantum-mechanical marginal distributions in phase space. (Published in J. Phys. A: 
Math. Gen. 36, 12275 (2003)) 



I. INTRODUCTION 



The concept of entanglement naturally appears in quantum mechanics when the superposition principle is applied 
to composite systems. In this sense, a multipartite system is entangled when their physical properties cannot be 
described through a tensor product of density operators associated to their different parts which constitute the 
whole system. An immediate consequence of this important effect has its origin in theory of quantum measurement 
refl: the entangled state of the multipartite system can reveal information on its constituent parts. However, this 
information is extremely sensitive to the dissipative coupling between the macroscopic meter and its environment. 
In fact, entangled states involving macroscopic meters are rapidly transformed into statistical mixtures ofproduct 
states and this fast relaxation process characterizes the decoherence According to Raimond et al. "The 

decoherence itself involves entanglement since the meter gets entangled with its environment. As the information leaks 
into the environment, the meter's state is obtained by tracing over the environment variables, leading to the final 
statistical mixture. This analysis is fully consistent with the Copenhagen description of a measurement" . Beyond 
these fundamental features, entangled states have potential applications for information processing and quantum 
computing (6j, iTj, jS, iSl] j quantum teleportation ilfll], dense coding 11], and quantum cryptographic schemes jl2l] . 

A feasible physical system to generate entangled states is given by the Jaynes-Cummings model (JCM) which 
describes the matter-field interaction 0, . It is typically realized in cavity QED experiments involving Rydberg 
atoms crossing superconducting cavities (one by one) in different frequency regimes and configurations, with relaxation 
rates small and well understood Sj. Recently, many authors have inve stig ated the two-mode and driven JCM in 
different contexts and predicted new interesting results [IllIlllllllllllMllllllllllMllllllll^. For instance, 
Solano et al. [2^ have proposed a method of generating multipartite entanglement through the interaction of a system 
of N two-level atoms in a cavity of high quality factor with a strong classical driving field. Following the authors, the 
main advantage of this external field in the system under consideration is the great flexibility in generating entangled 
states, since it provides freedom in choosing the detuning and strength of the field. On the other hand, Wildfeuer 
and Schiller zr\ have used the Schwinger's oscillator model to obtain a mathematical solution for the generation of 
entangled A^-photon states in the framework of the two-mode JCM. Here we develop a mathematical procedure which 
permits us to obtain compact solutions for atomic inversion and Wigner function in the framework of the driven 
JCM, considering any cavity and external fields. In particular, both solutions are expressed in the integral form with 
their integrands presenting a common term that describes the product of the Glauber-Sudarshan quasiprobability 



FIG. 1: Experimental apparatus used in the description of the JCM with an external quantum field. 



distributions |28| for each field, and a kernel responsible for the correlations. To illustrate our results we fix the cavity 
field in the even- and odd-coherent states E^li a-nd the driving field in the coherent state. Furthermore, we show how 
the detuning and amplitude of the driving field modify the entanglement in the tripartite system via Wigner function. 

The paper is organized as follows. In Section II we obtain the time-evolution operator and the matrix elements of 
the density operator for the driven JCM, with the cavity and external fields described in the diagonal representation 
of coherent states. Following, we fix the cavity field in the even- and odd-coherent states and the driving field in the 
coherent state to investigate, in Section III, the effects of amplitude of the driving field and detuning parameters upon 
the atomic inversion. In Section IV we derive a formal expression for the Wigner function associated with the cavity 
field which permits us to analyse how the entanglement is modified in the tripartite system. Moreover, we also obtain 
analytical expressions for the correct quantum-mechanical marginal distributions in phase space. Section V contains 
our summary and conclusions. Finally, Appendixes A and B contain the main steps to calculate the atomic inversion 
and Wigner function, respectively. 



II. ALGEBRAIC ASPECTS OF THE JCM WITH AN EXTERNAL QUANTUM FIELD 

In general, the driven JCM consists of a two-level atom interacting nonresonantly with a single-mode cavity field, 
and driven additionally by an external field through one open side of the cavity (the experimental scheme in the 
context of cavity QED is sketched in figure 1). Within the dipole and rotating-wave approximations, the dynamics of 
the atom-cavity system is governed by the Hamiltonian H = Hp + V, where 

Ho = hoj (a.^a + h%) + ^hujcr^ , (1) 
V = -M CTz -f ?iKa (a^<T_ -f acr_|_) -I- fi-Kb (b^(T_ -|- bcr_|_) . (2) 

Here, to is the cavity field frequency (we assume the resonance condition between the cavity and driving fields), ujq is the 
atomic transition frequency, S = luq — lu is the detuning frequency, and Ka(b) is the coupling constant between the atom 
and the cavity (external) field. The atomic spin-flip operators (t± and cr^ are defined as cr_|_ — |e)(5|, (T_ = |(7)(e|, and 
CTz = |e)(e| — \g){g\ {\g) and |e) correspond to ground and excited states of the atom), with the following commutation 
relations: [crz,cr±] = zL2(t± and [cr_|_,cr_] = a^- Furthermore, a (a^) and b (b^) are the annihilation (creation) 
operators of the single- mode cavity and external fields, respectively. It is important to mention that the quantum 
nature of the fields used in many proposed schemes for quantum information processing present serious consequences 
in large scale quantum computations, since the uncertainty principle and the possibility of becoming entangled with 
the physical qubits represent possible limitations on quantum computing |3(1 l3lL 133. ISa. l34l |. In this sense, van Enk 
and Kimble 31] have considered the interaction of atomic qubits laser fields and quantify atom- field entanglement 
in various situations of interest where, in particular, they found that the entanglement decreases with the mean 
number of photons (n) in a laser beam as E oc (n)~^log2(n) for (n) » 1. Pursuing this hne, Gea-Banacloche [s^l 
has investigated the quantum nature of the laser fields used in the manipulation of quantum information, focusing 
especially on phase errors and their effects on error-correction schemes (for more details, see [sslls^ l. 

Now, let us define the quasi- mode operators A = eaa -I- ebb and B = e^^SL — ejo (where ea(b) = '*a(b)/'*cff i Sind 
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^cff 



is an effective coupling constant), which satisfy the commutation relations 



[A, At] 


= 1, 


[Na,A] 




[NA,At] = At 


[B,Bt] 


= 1, 


[Nb,B] 


= B, 


[NB,Bt] =Bt 


[A,B] 


= 0, 


[A,Bt] 


= 0, 


[Na,Nb] = 0, 



being Na = AtA (Nb = BtB) the number operator related to the quasi-mode operator A (B). Introducing the 
number-sum operator S = Na + Nb and the number-difference D = Na — Nb, we verify that: 

(i) S = Ha -l- nb is a conserved quantity (na = ata and nb = btb are the photon- number operators of the cavity and 
external fields); 

(ii) the operator D can be written in terms of the generators {K+, K_, Kq} of the SU(2) Lie algebra, 

B = 2{el- el) Kq + 2e,eb (K_ + K+) , (3) 

where K_ = abt, K+ = atb, and Kq 5 (ata- btb), with [K_,K+] = -2Ko and [Ko,K±] = ±K±; 

(iii) the commutation relation between the operators S and D is null, i.e., [S,D] — 0; and consequently, 

(iv) the Hamiltonian H simplifies to 

Hq = tuj S + —huj (Tz , (4) 
V = iM(T^ + ftKcff(Atcr_ +Acr+) , (5) 

with [Ho, V] = 0. This fact leads us to obtain the Hamiltonian Hint = V in the interaction picture, which describes 
the well-known nonresonant JCM Hamiltonian for an atom interacting with the quasi-mode A, and whose coupling 
constant is given by Kcff- Thus, the unitary time-evolution operator is the usual nonresonant JCM time-evolution 
operator. 

If one considers the time-evolution operator I4{t) — exp (~iVt/h) of the atom-cavity system written in the atomic 
basis, the elements l^ij(i) of the 2x2 matrix can be expressed as 



ijof 

m 



Uii{t) = cosUV/^A . (6) 

^ ^ ^ V Pa 

, s sin (t^fBT ) , , 

U^2{t) = -iACeff }^ ' A , (7) 

. t sin(tV^' 

-IKoff AT Z— 



Z^2i(t) = -iKcffAT ^7=^, (8) 

vPa 

U,,[i) = cos(t7^)+i| ""(^Xg^^ , (9) 

where (^a = '^off-'^A + ^-nd /^a ~ Va + '^cff-'-- This result permits us to determine the density operator 

p(i) = Vl{t)p{0)U'' (t), being p(0) the density operator of the system at time t = 0. For convenience in the calculations, 
we assume the atom is initially in the excited state and the cavity and external fields are in the diagonal representation 
of coherent states, i.e., p(0) = Pat(O) ® Pab(O) '^ith Pat(O) = \^){^\ ^"^^ 

P.bm = PM®Pbm= // "% Pa(aa)Pb(Qb)|«a,ab)(aa,ab| , (10) 



where P{a) represents the Glauber-Sudarshan quasiprobability distribution for each field, and jaa, ah) = l^a) ® l^b). 
Consequently, the matrix elements pjj (t) can be calculated through the expressions 

PiiW - Unit) p,h{0)uUt) , (11) 
p,,{t) - lYn(t)Pab(0)WLW, (12) 
P2i(i) - lY2i(t)Pab(0)Wli(t), (13) 

P22W = l^2l(i)Pab(0)WkW- (14) 
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Note that p{t) describes the exact sohition of the Schrodinger equation in the interaction picture with the nonresonant 
driven- JCM Hamiltonian. Using this solution we can estabhsh analytical expressions for the time evolution of various 
functions characterizing the quantum state of the cavity field, such as the atomic inversion, the moments associated to 
the photon-number operator, the Mandel's Q parameter, the photon-number distribution and its respective entropy, 
the variances of the quadrature components, and the Wigncr function. For simplicity, the initial state of the driving 
field will be fix in the coherent state throughout this paper (Pb(0) — \0)(0\), and the initial state of the cavity field 
will assume two different possibilities: the even- and odd-coherent states [23 ■ With respect to Glauber-Sudarshan 

quasiprobability distribution, these considerations are equivalent to P^'^\ai,) = 7r(5^^^(ab — P) and 



4 cosh (lap) 

Trexp (|aaP) 
4 sinh(|ap) 



(5(2) (a^ -a) + 6('^ (a, + a) + 2 cosh (a^-a*^) J^^) J 
(5(2) (a^ -a) + 5(2) (a, + a)-2 cosh L-^-a*-^) S^'^ {a,) 



being S'^^^z) the two-dimensional delta function. According to Glauber " If the singularities of P{a) are of types 
stronger than those of delta function, e.g., derivatives of delta function, the field represented will have no classical 
analog" . Thus, in the next sections we will investigate the influence of the amplitude of the driving field and detuning 
parameters on the nonclassical effects of the cavity field where, in particular, the atomic inversion and the Wigner 
function should be emphasized. 



III. ATOMIC INVERSION 



The atomic inversion X{t) = Tr [p{t)(Tz\ is a quantity of central interest in this section since it is easily accessible 
in experiments 31]. For the atom-cavity system described in the previous section, this function can be written in an 
integral form as (see appendix A for calculational details) 



I{t) 



-Q!b 



-Pa(aa)-Pb(Q!b) ^iaa,,ah;t) 



(15) 



where 



S(aa, ttb; i) = 1 - 2 exp ^- |ea 



n=0 



n 



\Gn{tW 



Here, the function Gn{t) = — i(rin/A„) sin(A„t/2) is responsible for the time evolution of the atomic inversion, being 
— S'^ + ~ ^Kcffy/n + 1 the effective Rabi frequency. Note that the Eq. (|15|l can be obtained for any 

states of the cavity and external electromagnetic fields. For instance, if one considers the both cavity and external 
fields in the coherent states, the atomic inversion coincides with S(a,/3;t). This situation was investigated by Dutra 
et al. for the atomic excitation probability Pdt) = 5 [1(0 + 1] and S = (resonance condition), where the authors 
have shown that Peit) is connected to Wigner characteristic function of the cavity field since the conditions Ka 3> Kb, 
Kht <C 1, and \P\ 3> (Ka/Kb)tai (intense driving field) are satisfied. On the other hand, if one considers the cavity 
and external fields in the thermal and coherent states, respectively, the atomic inversion is given by 



Ithit) = 1- 



exp 



1 



E 

n=0 



1 



Inil 



\Gn{t)\' 



(16) 



In this expression, n is the mean number of thermal photons at time t — 0, and L„{z) corresponds to a Laguerre 
polynomial. Furthermore, the parameter ea(b) represents a scale factor for fij||/3|). It is important mentioning that 
Eq. H16|l corroborates the numerical investigations realized by Li and Gao 2lj for the thermal states, and this fact 
leads us to proceed with the study of atomic inversion for the even- and odd-coherent states. 



Let us consider the Glauber-Sudarshan quasiprobability distributions P^'^\aa.) and Pt"'{aa,) for the even- and odd- 
coherent states into the Eq. lfTB)l . whose integrals in the complex a^- and ab-planes can be evaluated without technical 
difficulties. In both situations, the atomic inversion is expressed in the compact form 



1 



t:\a,l3)\Gn{tt 



n=0 



2^;?(°)(a,/?)|G„(i)|^ 



(17) 



(18) 



n=0 
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FIG. 2: Time evolution of the atomic inversion Jc{t) of the atom initially prepared in the excited state interacting with cavity 
and external fields in the even-coherent and coherent states, respectively. These pictures correspond to (a,c) 5 — (resonant) 
and (b,d) S — Gkch (nonresonant) for ea = S/^/TO, eb = l/vTO, and |a| = 1 fixed, where two different values of amplitude of 
the driving field were considered: (a,b) j/3j = 2 and (c,d) |/3j — 20. 



with ^n\a, (3) and ^n' {a, P) given by 



exp (|c 



4cosh(|aP) 
+ 2exp(-2|a|2) Re 
exp (jap 



exp {- lea" + fb/^r 



2\ |eaa + eb/3| 



h exp y - \e^a - eb/3| j 



\e^a - eb/3|^ 
,1 



exp [(eaa + eb/3) [e^a - eb/3)*] 



(eatt + eb/3) (caa - eb/3) 



71! 

*1 ^^ 



4sinh(|Q;p) I 
-2exp(-2|anRe 



exp |eaa + eb/3|' 



2\ |eaa + eb/3| 



2n 



exp I eaa - eb/3|' 



2\ |eaa - eb/3|' 



exp [(eaa + eb/3) (eaa - eb/3)' 



[- (eaa + eb/3) (eaa - ehP)*] 



Fig. 2 shows the plots of Ic{t) versus Kcst when the atom-cavity system is resonant (a,c) (5 = and nonresonant (b,d) 
5 — 6Keff for ea — 3/vTO, eb = 1/a/10, and |a| = 1 fixed, with two different values of amplitude of the driving field: 
(a,b) |/3| = 2 and (c,d) |/3| — 20. Since the atom was initially prepared in the excited state, the value of the atomic 
inversion at the time origin is equal to one in all situations. In Fig. 2(a), we can perceive that Xo(i) behaves in a fairly 
irregular manner and the revivals are not well defined (in particular, the revivals are considered as a manifestation 
of the quantum nature of the electromagnetic field inside the cavity); while in Fig. 2(c), the collapses and revivals 
appear when the driving field is strong. Now, if one analyses the Figs. 2(b) and (d) we conclude that the collapses 
and revivals can be controlled by the detuning between the cavity (external) field and the atomic transition (in 
particular, the revivals have a regular structure and small amplitude). Similarly, Fig. 3 shows the plots oiXo{t) versus 
Kcfft considering the same parameter set used in the previous figure, where we verify that: (i) different structures of 
collapses and revivals are present, and (ii) the effects of the parameters |/3| and 6 on the atomic inversion Io{t) are 
completely analogous to the even-coherent states. Gora and Jedrzejek [s^ have shown that in the usual JCM with the 
cavity field prepared initially in a coherent state with a small mean number of photons (i.e., (n)c ~ 2 at time t = 0), 
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FIG. 3: Plots of Io{t) versus Kcfft G [0,200] for (a,c) S = (resonant) and (b,d) S — 6k,cs (nonresonant), with ta = 3/vTO, 
Eb ~ l/vTO, and laj — 1 fixed. In both situations were considered different values of amplitude of the driving field, i.e., (a,b) 
i/3i = 2 and (c,d) \/3\ = 20. 



the atomic inversion displays distinct collapses and revivals provided the atom and the field are slightly detuned, 
and the long-time behaviour of the model presents superstructures such as fractional revivals and superrevivals. In 
this sense, the Figs. 2(b) and 3(b) present a short-period behaviour with analogous superstructures and this fact 
is associated to the small mean number of photons used for both the cavity and external fields ((na(0))o ~ 0.762 
and (na(0))o ~ 1.313, with (nb(0))c ~ 4 fixed), since the detuning is large as compared to the effective coupling 
constant (e.g., S/2kcs = 3). Moreover, these superstructures disappear when we consider (nb(0))c ~ 400 in Figs. 
2(d) and 3(d). Summarizing, the amplitude of the driving field and the detuning parameter have a strong influence 
on the structures of collapses and revivals in the driven JCM, and this fact leads us to investigate its effects on the 
nonclassical properties of the cavity field via Wigner function. 



IV. WIGNER FUNCTION 

In many recent textbooks on quantum optics [35^ , the Wigner function is generally defined in terms of an auxiliary 
function (also denominated as Wigner characteristic function) which describes the symmetric ordering of creation 
and annihilation operators of the electromagnetic field, i.e., xiO = Tr[pD(^)] with D(^) = exp (^a^^ — ^*a) being the 
displacement operator. The connection between both functions is established by means of a two-dimensional Fourier 
transform as follows: 

Wh) = J^^^P - 7*0 xiO • (19) 

Thus, if one considers the cavity field in the framework of the driven JCM, its Wigner characteristic function can be 
defined in a similar form to atomic inversion, 

x(C;i) = JJ ^^^^f^Pa(aa)Pb(ab)K4(aa,ab;t) , (20) 
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FIG. 4: Plotsof VKi°'(7;t) versus p G [—7, 7] and q G [—10, 4] for the atom-cavity system with two different values of detuning: 
(a,b) 5 = (resonant) and (c,d) S — lOKcff (nonresonant), where the parameters |q| — 1 ({na)e ~ 0.762) and K^gt = 100 were 
fixed in the present simulation. In both situations, the condition Ka{b) = «^ was established and the values of amplitude of the 
driving field (a,c) |/3| = 2 ({nb)c = 4) and (b,d) 1/3| = 5 ((nb)c = 25) considered. 



with K^(Q!a, ab; t) given by 

K^{aa,ah;t) = (aa, ab|Wli(t)Da(C)i^ii(i)|aa, ^b) + (aa, ab|Z^2i(*)Da(C)^2i(i)l"a, ab) • 

Here, the displacement operator Da(C) is associated with the cavity field. Now, substituting x(C;^) ™to Eq. (|19|) . the 
expression for the Wigner function is promptly obtained, 

W,{r,t)^ ^^^^^f^Pa(«a)Pb(«b)K^(aa,ab;i) , (21) 
where the label 7 corresponds to representation in the complex phase-space and 

K^(aa, ab; = / ^ exp - 7*0 (aa, ^b; ■ (22) 

The functions Kj(aa,Q;b;^) and K-y(aa, ctb; i) were derived with details in the appendix B. In particular, when t — 
the function Kj{aa, ctb; 0) = 2 exp(— 2|7 — ckaP) does not depend of variables associated with the external field and 
this fact leads us to write the initial Wigner function as 

M^a(7; 0) = 2 / ^ exp(-2|7 - a,\^)PM . 
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FIG. 5: The Wigner function Wi (fji) is plotted assuming the same set of parameters estabhshed in the previous figure for 
the detuning frequency and ampUtude of the driving field, with |a| — 1 ({na)o ~ 1.313) and Kcst = 100 fixed. Note that the 
entanglement is maximum when ^ = (resonant regime), and minimum for S = lOACeff (nonresonant regime). 



This expression represents a Gaussian smoothing process of the integrand Pa(aa) such that Wa(7; 0) is a well-defined 
function in the phase space p = \/2Ini(7) and q = •\/2Re(7). On the other hand, for t > the function K^(aa, ah', t) 
is responsible for the entanglement between the cavity and external fields (here represented by the Glauber-Sudarshan 
quasiprobability distributions Pa(Q^a) and Pb(ab), respectively) since the complex variables a^, and at are completely 
correlated. Furthermore, it is important mentioning that xi^it) and VFa(7;t) can be evaluated for any states of the 
cavity and external fields (similar condition was established for atomic inversion) without restrictions on the different 
interaction times, and the ex pres sions obtained analytically from this procedure generalize the results previously 
discussed in the hterature [1^ [13 ■ 

For instance, let us consider the Glauber-Sudarshan quasiprobability distributions for even- and odd-coherent states 
into the Eq. H21(l . After the integrations in the complex a^- and ab-planes, we get 

^) = [K7("> t) + K^(-«> /3; t) + exp(-2|aniK^(a, /3; t)] (23) 

4cosh(|ap) 



and 



where 



Wt\rt) ^ /'^'^T^'f |l [K-,(a,/?;0+K^(-«,/?;^)-cxp(-2|a|^)IK,(a,/3;t)] , (24) 
4smh(|Q;|-^) 



iy{a, P; t) 



( d d \ 

2exp(|Qn cosh ( -"*^ ) [exp(|aanK^(aa,/3;0] 



(25) 
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Note that at time t = 0, the function exp(— 2|ap)K^(a, /3; 0) = 4exp(— 2|7p) cos[4Ini(7Q;*)] leads us to recover 
weU-known expressions in the hterature p9i |: 

Wt\r,0) = ^^^exp(-2|7n {exp(-2|ancosh[4Re(7a*)] +cos[4Ini(7a*)]} 



and 



Wt\T,0) = ^^^M^exp(-2|7n {exp(-2|an cosh[4Re(7«^ 



— cos [4Im(7a*)]} 



Figs. 4 and 5 show the three-dimensional plots of Wa'^\'j;t) and Wa°\j-t) versus p = \/2Im(7) and q = y2Re(7), 
respectively, for the atom-cavity system resonant (a,b) S = and nonresonant (c,d) 5 — lOKcff- In both simulations, 
we consider |q;| = 1 and two different values of amplitude of the driving field: (a,c) |/3| = 2 and (b,d) \(3\ = 5. 
Furthermore, we also fix the parameter K^gt — 100 which permits us to obtain a partial view of entanglement in the 
tripartite system. A first analysis of these pictures shows, via Wigner function, that the entanglement is sensitive to 
variations of the experimental parameters |/3| and 6 (this fact corroborates the previous results obtained for atomic 
inversion); being the detuning parameter responsible for the entanglement degree between the components involved 
in the system, since both the driving and cavity fields arc in resonance. In this sense, although the Figs. 4(c)-(d) and 
5(c)-(d) have similar structures, there are subtle differences between them: Wa°H^]t) assumes negative values due to 
initial sub-Poissonian photon statistics of the cavity field; while wi''^ (7; t) is strictly positive, since the even-coherent 
state has super- Poissonian photon statistics for any initial value of (na)e (see Ref. ,29i| for more details). On the 
other hand, the increase of in Figs. 4(b,d) and 5(b,d) shows an interesting effect on the Wigner functions: the 
interference patterns between the states of the driving and cavity fields turn to be more pronounced, and this effect 
modifies the shapes of wi°'' (7; t) and wi""* (7; t). Similar analysis can be also applied if one considers both the external 
and cavity fields in the coherent states (see appendix B). 

To conclude this section, we will determine the marginal probability distribution functions {ipniq; t)\'^ and \(pii{p; i)P 
through the direct integration of Eq. H21|l over the variables p or q, i.e.. 



IV'a(g;i)P = 
lva(p;i)l' = 



^'a(aa)^b(ab)Ug(Q!a,ab;i) , 



^'a(cka)Pb(ab) Vp(aa, ab; t) , 



(26) 
(27) 



with 



U,(aa,ab;i) 



dp 

7^ 



K^(aa, ah] t) 



and 



Vp(aa,ab;i) 



dq 

7^ 



For this purpose, let us initially introduce the complex function 

{m.m } 



H('"''"')(aa,ab) = J2 (2fc)!!4""'^(0)4'"'"'^(0) 



k=Q 



£a"a + Eb^b 



a/2 



eaeb(ebaa 
V2 



CaOfb) 



where Hn{z) is the Hermite polynomial, fa(b) = ea(b) (ea(b)aa 
to'. In addition, we define the auxiliary functions 



eb(a)Ctb)i a-iid {m, m'} stands for the minor of m and 



Y(™'"')(aa,ab;t) = 7i(:"'"''(aa,ab)F„(t)F*,(i) 



CaEb ebQ!a 



+ ga«b .^(m+l,m/ + l) ^^-j Gm{t)G*^,{t) 

2 Ea^a + CbCkb ^ {m + 1)(to' + 1) 



and 



4™''"')(aa,«b) = V2 



exp 



V2 



ebttbl 



[v/2eb (ebCKa + Caab)]"" [\/2ea (ea^a -I- Eb^b)*]' 

(2m)!! (2to')!! 
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which permit us to express the integrands Ug(aa, ctb! t) and Yp(ag,, ctb! t) in compact forms as follows: 



oc 



Ug(aa,ab;i) 



(28) 



m,rn'—0 



oo 



Vp(aa,ab;i) 



iaa,-iab) 



itta, -iab; t) . 



(29) 



m,rn'—0 



Consequently, the connection between Eqs. H28() and (|29|l can be promptly established through the mathematical 
relations \Jg(a^, ab;t) = Yg{ia^,iab;t) and Vp(aa, ab; *) = Up(-iaa, -iftb; *)■ 

In analogy to Wigner function, the marginal probability distribution functions do not depend on the driving field 
at time i = 0, since their expressions are reduced to 



Now, if one considers both the external and cavity fields in the coherent states, we get |V'a('?;OP = Uq(a,/9;<) and 
|(/7a(p;0P — Vp(a,/?;t). In particular, this example shows that the marginal distributions represent an important 
additional tool in the qualitative study of entanglement, since the variables a and /3 are completely correlated. 



In this paper, we have applied the decomposition formula for SU(2) Lie algebra on the driven Jaynes-Cummings 
model in order to calculate, for instance, the exact expressions for atomic inversion and Wigner function when the 
atom is initially prepared at the excited state. In fact, adopting the diagonal representation of coherent states, we 
have shown that these expressions can be written in the integral form, with their integrands presenting a commom 
term which describes the product of the Glauber-Sudarshan quasiprobability distribution functions for each field, and 
a kernel responsible for the entanglement. It is important mentioning that the mathematical procedure developed 
here does not present any restrictions on the states of the cavity and driving electromagnetic fields. Following, to 
illustrate these results we have fixed the driving field in the coherent state and assumed two different possibilities 
for the cavity field (i.e., the even- and odd-coherent states). In this way, we have verified that the amplitude of the 
external field and detuning parameter (i) perform a strong influence on the structures of collapses and revivals in the 
atomic inversion, (ii) control the entanglement degree in the tripartite system; and consequently, (iii) modify the shape 
of Wa(7; t) since the interference patterns between the states of the driving and cavity fields turn to be more evident 
through the Wigner function. In addition, the formalism employed in the calculation of atomic inversion and Wigner 
function open new possibilities of future investigations in similar physical systems (e.g., see Refs. H^H^); or in the 
study of dissipative composite systems, where the decoherence effect has a central role in the quantum information 
processing. These considerations are under current research and will be published elsewhere. Summarizing, the work 
reported here is clearly the product of considerable effort and represents an original contribution to the wider field of 
entangled-state engineering with emphasis on quantum computation and related topics. 
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APPENDIX A: THE INTEGRAL FORM OF THE ATOMIC INVERSION 



With the help of the definition estabhshed in Sec. Ill for atomic inversion and the cychc invariance property of the 
trace operation, we get 



m - Trab{pab(0) [wIi(i)WiiW-Z^JiWW2i(t)]} 

J2 sin2 (ty^ ) 



Trab <^ Pab(O) 



2 /3, 



(Al) 



Employing the diagonal representation of Pab(O) the coherent states basis into the second equality of Eq. (|A1|I . the 
integral form of the atomic inversion can be promptly obtained, i.e.. 



where 



At)^ II ^^^4^^aMPb(ab)S(aa,ab;i) 



j2 sin^ it\/f^A ) 
E{as,,ah;t) = (eta, abl cos(2iy^)|Q;a, ab) + -^r (aa,ab| |aa,ab) 



(A2) 



(A3) 



However, the effectiveness of the integral form (|A2|) is connected with the determination of an analytical expression 
for Eq. (K^ . 

To calculate the function S(Q;a, a\,; t), firstly we expand the operators cos(2t-\/ /J^ ) and sm'^{ty/ (3j^ ) / in a power 
series as follows: 



k=0 



x=0 



and 



sin2(ty^) _ 1 — ^-1^*= 



^[2(^ + 1)1!^'^^''^ dx^ 



x=Q 



Secondly, we apply the antinormal-order decomposition formula for SU(2) Lie algebra on the operator e^^, which 
leads us to obtain [sllsgllioll 



^xD ^ gB+K_ gS+SoK+ g(lnBo)Ko 



where 



2eaeb sinh x 



cosh.T + (e^ — e'^) sinh.T 



and Bq = [cosh a; + (e^ — e^) sinh.T] ' 



After lengthy calculations, the analytical expressions for the mean values 

oo 

(tta, Ckbl C0s(2ty^)|aa, Ckb) = CXp lea^a + Ebttbl^) ^ 



Ea^a + Cbttb 



l2n 



n 



■ cos(t A,i 



(A4) 



and 



,sin2(tv/^), 
/3a 



l^a, ah 



) = exp leaQa + fbabl^) ^ 



2\ Ifa^a + Cbabl^" sin^(tA„/2) 



nl (A„/2)2 



(A5) 



are determined, with ^ S"^ + flf^ and fin = + 1 the effective Rabi frequency. Now, substituting these 

results into Eq. IjASp . we obtain 



oo 

S(aa, Q!b; = 1 - 2 exp ^- [eatta + ebCKbl^) ^ 



(A6) 



n=0 



where G„(i) = — i(ri„/A„) sin(A„t/2). Consequently, with the determination of the analytical expression for 
S(aa,Q;b;i)i the effectiveness of the integral form (|A2p is guaranteed. 
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APPENDIX B: CALCULATIONAL DETAILS OF THE WIGNER FUNCTION 

Initially, we will derive the means values 

(«a,ab|wLWDa(e)Z^ii(<)|aa,ab) - JJ '''y^'^ (aa,ab|Z^Ii(t)Da(OI/3a,/3b)(/3a,/3b|^^ii(OI«a,ab) (Bl) 



and 



by means of integrations in the complex variables /3a and /3b- Thus, let us substitute into Eqs. (|Bip and (|B2|) the 
auxiliary mean values 



(aa,ab|Z^^i(i)Da(e)Z^2iW|aa,«b) = ^ ^^^f^ (aa, ab|Z^^i WDa(OI/3a, /?b)(/?a, /3b|i^21 Wl^a, ab) , (B2) 



(aa,ab|i^IiWDa(e)|/3a,/3b) = exp - (^/3* - ^/^a) ((/3a + e,/3b|l^iiWI«a,ab))* 
(aa,ab|wLWDa(OI/3a,/3b) = exp i (^/S: - T/^a) ((/?a + C,/3b|Z^2iW|aa,ab»* 

oo 

(/3a,/3b|Wii(i)|aa,ab) - ^F„(i)A (aa,ab,^a,/^b) , 



m=0 



(/3a,/3b|W2l(i)|aa,ab) = (Ca/^a + eb/3b)* 



Vm + 1 



Am(aa,ab,/^a,/3b) 



with 



A„i(aa,ab,/3a, ^b) = exp 



"2 + l^bP + |/3aP + l/SbP) + (Cb^a " eaOb) (cb/^a " Ca/^b)* 



[(eatta + Cbab) (ea/3a + eb/3b)*] ' 



and Fm{t) — cos(A„ii/2) — i((5/Am) sin(A„ii/2) (the function Gm{t) was previously defined in appendix A). Then, 
carrying out the integrations in the variables /3a and /3b, we get 



and 



,ah\ll\i{t)'Di,{^)Un{t)\aa,ah) ^ ^ T^™'" ^ (cia, ab)4'"''" ^Oa, "b; 



m.m'—O 



(aa,ab|WkWDa(0i^2i(i)l«a,ab)= E T(™''")(aa,ab)jf''"^( 



(B3) 



(B4) 



where 



T> '(aa,ab) = exp 



1^12 ^ ^ 

lea^a + Ebttbl^ + Eb (eb^a - ea^b)* C ^ ^a (CaCka " eb^b) T 



[ea(eaaa + ebab)*^] 



m'-m lea^a + Ebttbl 



4™''"'^(aa,ab;i) = lL"''-") 
(aa,ab;i) = 



-j(m,m ) 



eb«a + Eattb 1^12 

fafb 1^1 

EaOa + EbCKb 



m'! 



m+1 {m'-m) 

m' + 1 ™+i 



eaCb ■ ?| 

CaQa + Ebttb 



Consequently, the function Kj(aa, Q!b; t) which appears in the integrand of t) can be determined as follows: 



K^(aa,ab;i)= ^ T^™'" ^ (cka, ab)r^'"'™ ^Oa, ab; i) , 



(m.m') / 



(B5) 



m.m'— 
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FIG. 6: Plots of Wa, {'y;t) = K^(ci,/3;t) versus p G [—7,7] and q G [—10,4] for the atom-cavity system resonant (a,b) 5 = 
(maximum entanglement) and nonresonant (c,d) S = lOfieff (minimum entanglement), with jaj = 1 and ACcfft = 100 fixed. We 
also have considered two different values of amplitude of the driving field: (a,c) = 2 and (b,d) |/3| = 5, where the condition 
'^a(b) = K was established in both situations. 



, . -p,(m,Tn')/ ,N -rfm.m')/ ,n -,(m,m')/ ,n 

being '(aa,ab;i) = ' (aa, ^b; + '{aa,ah;t). 

An immediate application of this result is the calculation of K^(Q!a, a^; t) since both functions are connected by a 
two-dimensional Fourier transform. Now, substituting IjBSII into Eq. (|22|) and integrating in the complex variable ^, 
we obtain as result the analytical expression 



where 



and 



K^(aa, ah] t) = 



m^rn'—Q 



C(™''"')(aa,ab)M('"'™')(aa,ab;i) , 



C^^'^'^aa, ab) = 2exp(-|eaaa + ebab|^-27a7^) [2ea (ea^a + eb^b)* 7a 

[(^a ~ ^b) (^aOa " eb^b) (Catta + Ebab)*]" 



m'l 



M(™'"')(aa,ab;i) = L^-""'^ 



4ea£b £bQa + £a«b 
^a - ^aCta - Cb^b 



7a7b 



(B6) 



to' -I- 1 ^ 



Eatta + Eb^b 



(m' — m) 
m+1 



4eaeb £bQ!a + eaCtb 

^a - ^b '^atta " eb^b 



7a7b 
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with 7a(b) = 7 — ea(b) (ea(b)aa — eb(a)«b) ■ In particular, when Ka(b) = k the expression for K-y(aa, ab; t) can be written 
in the simphfied form 



K^(aa,ab;i)- ^ 0(,™)(aa, ab) [o(,"')(aa, ab)] * 7^^™^'"'Haa, ^b; , 



(B7) 



m.m'— 



where 



0('")(aa,«b) 
7^(™''"')(aa,ab;^) 



= x/2 



exp 



4(1 



tta + ttbl +|27-(Q;a 



abj 



{(aa + Cth) [27 - (tta - "b)]*}' 



2"m! 



ab) 



V(m + l)(m' + l) 



This solution is equivalent to consider that the interaction between atom and cavity (external) field has the same 
strength. Note that (|B6(I represents an important step in the process of investigation of the effects due the amplitude 
of the driving field and the detuning parameter on the nonclassical properties of the cavity field via Wigner function. 

For instance, when the external and cavity fields were described by coherent states, the Wigner function coincides 
with Kj{a,P;t) and for t = 0, we obtain the initial Wigner function wt\-/;0) = 2exp(-2| 7 — a\ ). Fig. 6 shows 

the three-dimensional plots of Wa'^\^;t) versus p — \/2lm{'y) and q = \/2Ke{'-/) considering the atom-cavity system 
resonant (a,b) 6 — and nonresonant (c,d) 6 — IOkcS for \a\ — 1 ((na)c — 1) and Kcfft — 100 fixed, with two different 
values of amplitude of the driving field: (a,c) |/?| = 2 ((nb)c — 4) and (b,d) \(3\ = 5 ((nb)c = 25). The condition 
'*a(b) = ^ "ws-s established in both situations, and the infinite sums present in (|B7|I were substituted by finite sums as 
follows: 



771 — 



E E 

■m— m'— m+1 



O^") (a, (3) [0(™') (a, (3)] * 7^(™'™') (a, /3; t) 



where £ is the maximum value which does guarantee the convergence of this expression (we have fixed € = 50 in 
the numerical investigations). Since the time evolution of composite systems leads us to the essential concept of 
entanglement (Ej, the Figs. 6(a)-(d) reflect the effects of the driving field on the different forms of entanglement in 
the tripartite system for a specific value of Kost (maximum entanglement when S — 0, and minimum entanglement for 
S — lOKcff). For a global view of entanglement of the system under consideration, different values of Kcgt and (|/3|, 6) 
are necessary. Here, we give only a partial view of this important effect. 
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